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1. Introduction 

Over the past decade, model atmospheres and synthetic spectra for late-type stars 
have improved hand-in-hand with higher quality opacities. In 1994, quality lists of tran- 
sitions of the water vapor molecule based on ab initio molecular calculations became 
available (Miller et al. 1994; Schryber, Miller, & Tennyson 1995; Jorgensen, Jensen, & 
Sorensen 1994) which allowed the computation of the first Direct Opacity Sampling (here- 
after dOS) model atmospheres for late-type dwarfs (Allard et al. 1994) and brown dwarfs 
(Allard et al. 1996), later to become the NextGen models described in (Hauschildt, Allard, 
& Baron 1999). Showing a more physical description of their main opacities, the NextGen 
dOS model atmospheres promised a better description of the Spectral Energy Distribu- 
tion (hereafter SED) of cool stars. And this appeared to be verified for the infrared SED of 
M dwarfs Qones et al. 1996; Leggett et al. 1996; Allard et al. 1997). 

But despite these fundamental improvements, the NextGen models have failed to 
match adequately several of the optical (spectroscopic and photometric) properties of 
late type dwarfs and giants. In fact, the dOS models (Brett 1995; Allard et al. 1994, 1997; 
Hauschildt, Allard, & Baron 1999) provided a worse fit to the optical SED of lower main 
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sequence dwarfs than previous models based on simplijfied molecular opacities (AUard 
& Hauschildt 1995, hereafter, AH95). The models could not reproduce the tight relation 
formed by M dwarfs in the V-R versus R-I two-color diagram, indicating a systemati- 
cally shallower slope of the optical SED (defined by TiO absorption) then observed in 
these stars. A systematic flux excess in the spectral region sampled by the V bandpass 
(0.4 — 0.65/im) was noted as well in dwarfs as in red giants. Baraffe et al. (1998a) ob- 
served that this excess in the NextGen dwarf models translated into lower main sequence 
isochrones deviating progressively to the blue (by up to 1.0 magnitude!) in (Mi^,V-I) color- 
magnitude diagrams, for masses lower than about 0.5 M© (Teff < 3800 K). Baraffe et al. 
(1997a) examined a variety of globular clusters and showed that this departure of the 
models decreased in amplitude with decreasing metallicity. The problem seemed there- 
fore confirmed to be caused by a lack of opacity of an oxygen compound^. 

Only three independent models of the TiO molecule and corresponding lists of tran- 
sitions were available so far to the construction of model atmospheres. The first model 
was constructed over two decades ago by Collins (1975) and was restricted by the compu- 
tational limitations of the time. The Collins line list was intended to model the extended 
atmospheres of red giants and did not include high energy and otherwise weak transi- 
tions important by their number in the hotter environments of red dwarf atmospheres. It 
also neglected the red e system of TiO. Jorgensen (1994) extended Collins' work to TiO iso- 
topic transitions, included the e system from revised molecular rotational constants, and 
adopted the laboratory oscillator strengths of Davis, Littleton, & Phillips (1986). It is there- 
fore understandable that the resulting limited list of transitions caused shortcomings in 
the NextGen model atmospheres. The second TiO list was constructed by Kurucz (1993) 
and is used in his ATLAS9-12 atmospheres. The third model was constructed by Plez 
(1992) using also the Davis, Littleton, & Phillips (1986) oscillator strengths, and is used in 
his version of the MARCS atmosphere code. All three independent models yielded the 
visual flux excess in different proportions. Plez (1998) suggested that the missing opacity 
is due to missing TiO band systems in current lists, and added the TiO a-f system at 0.5 
/xm to his list. However, this helped him only partially to resolve the V-band flux excess 
problem. 

Recently, Langhoff (1997) constructed a new model of the TiO molecule and pub- 
lished new lifetimes and oscillator strengths that improved significantly upon the 1986 
values of Davis et al. (Valenti, Piskunov, & Johns-KruU 1998; Plez 1998). Schwenke (1998) 
has subsequently computed a corresponding list of transitions complete to the high ener- 



^ Hydride absorption bands only get stronger relative to the continuum with decreasing metallicity over 
the range covered by the globular clusters studied in Baraffe et al. (1997a) in the optical spectral range. 
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gies and therefore more suitable for general model atmosphere applications. In this paper 
we present the results of including this new TiO line list as well as the new AMES-H2O 
list in the construction of model atmospheres and synthetic spectra for late-type dwarfs 
and red giants. 



2. Model calculations 

We have calculated the models presented in this paper using version 10.3 of our gen- 
eral model atmosphere code PHOENIX. Details of the code and the general input physics 
are discussed in Hauschildt, AUard, & Baron (1999) and references cited therein. The 
models for M giants were calculated with the same setup, however, they employ spheri- 
cal geometry (including spherically symmetric radiative transfer). For giant models with 
low gravities (log(g) < 3.5), this can be an important effect for the correct calculation of 
the structure of the model atmosphere and the synthetic spectrum (Aufdenberg et al. 
1998, 1999). The main difference between the models presented in Hauschildt, AUard, & 
Baron (1999) and the models presented here is the use of the new AMES line lists for H2O 
(Partridge & Schwenke 1997) and TiO (Schwenke 1998), but we have also adjusted the 
empirical oscillator strengths of VO and CaH absorption bands to respect their strength 
relative to TiO bands (note that VO and CaH absorption is still treated in the Just Over- 
lapping Line Approximation due to lack of adequate line data). Our combined molecular 
line list includes about 500 million molecular lines. These lines are treated with a dOS 
technique where each line has its individual Voigt (for strong lines) or Gauss (weak lines) 
line profile (in the standard OS method tables of precomputed opacities are used). They 
are selected for every model from the master line list to the beginning of each model it- 
eration to account for changes in the model structure, see Hauschildt, AUard, & Baron 
(1999) for details. This procedure selects about 215 million molecular lines for a typical 
giant model with T^n ~ 3000 K and about 130 million molecular lines for a dwarf model 
with the same effective temperature. Therefore, we generally use the parallelized version 
of phoenix (Hauschildt, Baron, & Allard 1997; Baron & Hauschildt 1998; Hauschildt & 
Baron 1999) to perform the calculation efficiently on parallel supercomputers. Details of 
the TiO and H2O lists are given in the next subsections. 



2.1. Water lines 

The effects of water lines on the M dwarfs SED was discussed in Allard et al. (1994). 
For the work presented here, we have replaced the UCL water vapor line-list (Miller et al. 
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1994; Schryber, Miller, & Tennyson 1995, hereafter: MT-H2O) used in Hauschildt, Allard, 

& Baron (1999) with the AMES water line-list (Partridge & Schwenke 1997, hereafter: 
AMES-H2O). This list includes about 307 million lines of water vapor. For the calcula- 
tions shown in this paper we have used Hj^O and neglected other, much less abundant, 
isotopes of this molecule. 

The water vapor opacity is governed by the completeness of the line list used, but also 
by the adopted atomization energy. The partition function of the molecule cancels out in 
the final absorption coefficient, after we have multiplied cross-sections by number densi- 
ties. But since water is an important chemical equilibrium specie, errors in the partition 
function can affect indirectly the model structure and spectra. The AH95 models were 
based on the Ludwig (1971) hot flames water cross-sections in the form of straight means, 
and used the JANAF partition function for water vapor (Irwin 1988). The NextGen mod- 
els where, on the other hand, computed with the MT-H2O line list and a partition func- 
tion computed from the MT-H2O levels. We note that the AMES-H2O partition function 
is practically identical to JANAF values, while the MT-H2O value is smaller than JANAF 
for temperatures above 3000K, possibly due to the energy levels missing in the MT-H2O 
data. We have therefore adopted for this and later work the JANAF partition function. 
We use an atomization energy of 9.5119 eV from Irwin (1988) for all models since AH95. 



2.2. TiO lines 

The main point of this paper is the comparison of the model structure and the syn- 
thetic spectra obtained by using the list of TiO lines from (Jorgensen 1994, hereafter: 
SCAN-TiO) and the new list of TiO lines from (Schwenke 1998, hereafter, AMES-TiO). 
The AMES-TiO list includes a total of about 172 million lines, about 44.6 million of these 
are for the most abundant isotope "^^TiO and about 32 million lines for each of the re- 
maining 4 isotopes beyond the completeness of the line-list, two more 
considerations affects the overall opacity produced by TiO, and explain systematic dif- 
ferences between model versions and by different authors: the atomization energy (Dq) 
determines the number density of TiO, and the TiO band oscillator strengths^ have been 
derived from sunspot observations (Davis, Littleton, & Phillips 1986, hereafter: DLP86), 
laboratory experiments (Hedgecock, Naulin, & Costes 1995, hereafter: HNC95), as well 
as from ab initio calculations (Langhoff 1997, hereafter L97). (Brett 1990, hereafter: B90) 



^ffi ~ fi/v"Mv'p", where the /'s are the oscillator strengths and the q's are the Franck-Condon factors of 
the transition v'u" 
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derived astrophysical fei values by jfitting the optical SEDs of red giants, using an atom- 
ization energy of 7.76 eV. He quoted that reducing this value by 0.3 eV would increase his 
fei by a factor 2.5. The most recent estimate of Dq for TiO is now 6.92 eV, which suggests 
that the B90 fei values are underestimated by as much as a factor 7! We summarize in table 
1 the various sources of oscillator strengths available for TiO. 

The models of AUard (1990) used fei values from B90 together with the straight-mean 
TiO opacities by Collins (1975) and Collins & Fay (1974), and assuming an atomization en- 
ergy of 6.87 eV. The first comparison of these models to the SED of M dwarfs (Kirkpatrick 
et al. 1993) revealed the inadequacy of this combination of parameters for TiO which pro- 
duced far too weak optical opacities. We have therefore, since the AH95 model series, 
employed the updated value of 6.92eV, together with the larger laboratory fei- values of 
DLP86. These two modifications combined to significantly increased the strength of TiO 
opacities in the models, bringing the AH95 and later the dOS NextGen models in im- 
proved agreement with the SED of M dwarfs. Any differences in the predictions of the 
AH95 and NextGen models are therefore purely due to the opacity technique (Straight- 
Mean versus dOS) and to the completeness of the line-list used. The incompleteness of the 
SCAN-TiO line-list allows photons to escape between absorption bands (see e.g. Valenti, 
Piskunov, & Johns-KruU 1998), and thus leads to systematically and increasingly (with 
higher Tgff) bluer optical colors (V-I) than observed (Baraffe et al. 1998b, 1997b). For the 
current models we therefore explore the use of the more complete AMES-TiO line-list, 
and the yet larger theoretical fei-values of L97. 



3. Results 

We have calculated a number of model atmospheres using either the SCAN-TiO or 
the AMES-TiO list of TiO lines and using either AMES-H2O or MT-H2O as source of the 
H2O lines. All the other input physics is the same for both sets of models. All mod- 
els have been fully converged with their respective set of parameters. Note that these 
models have been constructed for the purpose of this paper only and not to model indi- 
vidual stars and thus do not include dust formation and opacities which is important in 
atmosphere models with effective temperatures below about 2500 K, such models are pre- 
sented in a subsequent paper. In the following, we will discuss the results for the dwarf 
and giant models separately. The baseline for our comparisons are the NextGen models 
(Hauschildt, Allard, & Baron 1999) for the dwarfs and the NG-giant models (Hauschildt 
et al. 1999) for the giants. 
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3.1. Effects of different TiO line lists 

The models discussed in this section were all calculated using AMES-H2O to isolate 
the effects of different TiO line lists on the model spectra and structures. 



3.1.1. M dwarf models 

In Figs. 1 and 2 we show a comparison of model spectra calculated with AMES- 
TiO (full curves) and with SCAN-TiO (dotted curves) (both using our adopted fei-set, as 
quoted in table 1) for several effective temperatures. The gravity (log(g) = 5.0) and abun- 
dances (solar) were selected to be representative of M dwarfs in the solar neighborhood. 
In both figures, the resolution of the synthetic spectra was reduced by boxcar smoothing 
to 20A. At high effective temperatures, the two sets of models are nearly identical due to 
reduced importance of TiO absorption. At very low Teff the two line lists apparently agree 
very well since only the lowest levels of TiO remain populated. It is essentially between 
Teff ~ 2000 K and ^ 3500 K that the largest completeness and quality effects of the TiO line 
lists are seen. 

Fig. 3 indicates the location of each TiO band system for a 2900K model. From this 
it becomes clear that the addition of a-f transitions, which depress the continuum from 
0.4 to 0.5/xm, is one of the largest improvements brought by the AMES-TiO list to our 
models. We note that the entire optical regime from 0.4 to 0.75/im shows generally more 
opacity in the AMES-TiO models then using the Jorgensen (1994) (hereafter: J94) line list. 
The e bands at 0.82 to 0.88/im have a more precise shape in the AMES-TiO list, and come 
out stronger as well. This is a result of the completeness of the AMES-TiO list which also 
removes flux excess escaping between the troughs of the bands. We note however that 
some regions, such as in the 7 band near 0.78/um, show less opacity in the new models. 

The main effect of the new AMES-TiO on spectroscopic and photometric Tgff esti- 
mates will however be dominated by the change we make to the oscillator strengths. The 
L97 fei values being generally smaller than the DLP86 values adopted by J94, models of 
early-type M dwarfs using the new AMES-TiO setup should predict systematically lower 
effective temperatures then did prior models (NextGen, AH95, etc., see also Fig. 10 be- 
low). And beyond the enhanced completeness of the AMES-TiO list to high temperature 
transitions, the need for a cooler model should also contribute to making the TiO bands 
fit better a given star i.e. larger bands with less flux escaping from deeper, hotter atmo- 
spheric regions between them. 
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We could have opted to use the HNC95 laboratory fei-values as did Plez (1998), but 
since the L97 ab initio values agree quite well, we decided to keep these, except for the S 
and (fi band systems. The reason the oscillator strengths for the S and cp bands are less ac- 
curate is that it is very hard to get a good description of the b state, which is the upper state 
in both bands. For the S system, L97 derives an oscillator strength which is, as opposed 
to all other bands, twice as large as the DLP86 value. And the 6 and ip fepvalues cannot 
be corroborated by recent experimental values. Such a strong S band system would be 
difficult to bring in agreement with M dwarfs observations. Indeed, prior models have all 
shown a gradually increasing departure to the blue of the main sequence in My vs y — / 
diagrams (Baraffe et al. 1995, 1998b). Such departure is significantly improved using 
the new TiO list if one keeps a weak S band as indicate preliminary results of evolution 
models to be published separately (see also Fig. 12 below). We have therefore adopted to 
keep the DLP86 oscillator strength values for the reddest two TiO bands until new labo- 
ratory experiments can either confirm or infirm the L97 predictions. The summary of our 
adopted set of oscillator strengths for TiO is presented in Table 1. 



3.1.2. M giant models 

The results for the giant models are similar to the results for the dwarfs. Figures 4 
and 5 show synthetic spectra for 3 representative giant models with the indicated effective 
temperatures. The models have in common the parameters log(g) = 0.5, M = 5M0 and 
solar abundances. The differences between the AMES-TiO (full curves) and SCAN-TiO 
(dotted curves) models are somewhat larger for giants than for the dwarfs in the blue 
spectral region due to an increased sensitivity to the added a-f system opacities in the 
AMES-TiO line list. It is however somewhat less pronounced in the red spectral region 
where TiO bands are weaker in giants. The "spikes" that are apparent in the SCAN-TiO 
spectrum with T^h — 3000 K are absent in the AMES-TiO models. These spikes were one 
of our major problems in fitting observed spectra of giants. For larger Teff the differences 
between the spectra diminishes quickly as TiO becomes less important in the giants. This 
happens at lower effective temperatures compared to the dwarfs because of the lower 
pressures in giant atmosphere which results in smaller partial pressure of molecules as 
compared to dwarfs. 
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3.2.3. Model Structures 

A comparison of the model structures for both dwarf and giant model reveals only 
very small differences between structures calculated with AMES-TiO and SCAN-TiO. We 
plot the differences in electron temperature as well as the relative differences between the 
AMES-TiO and SCAN-TiO models for dwarfs and giants in Figs. 6 and 7, respectively. 
The changes are generally very small, only in the per-cent range for the gas pressures and 
about 10 K maximum difference between the electron temperatures for the dwarf models. 
For the giant models the differences are somewhat larger. The changes in the opacity 
averages are generally small but largest for the Rosseland mean opacity in the outer layers 
of the giant models. The temperatures are higher in AMES-TiO model for both the giant 
and the dwarf models, however, the gas pressures are lower in the AMES-TiO dwarf 
model but higher in the AMES-TiO giant model. Overall the changes are modestly small, 
indicating that the detailed effects of the TiO line lists do not have a large effect on the 
model structure itself. 



3.2. Effects of different water vapor line lists 

In figure 8 we show the effects of different water line lists on the synthetic spectra for 
M dwarfs. All models shown in the graph otherwise use the same line-lists (AMES-TiO 
was used for the TiO lines). In overall, we can see that the changes in the water vapor 
line lists are of larger amplitude than the changes in the completeness of the TiO line 
list. The models calculated with AMES-H2O show a totally different shape of the 1.4 /xm 
band, both weaker and wider then predicted by the MT-H2O model. The completeness of 
the new line list to high temperatures helps block more flux escaping from deeper, hotter 
layers of the models around 1.6 /xm, and 2.2 /xm. This promises a much better description 
of observations in general. 

We also find important changes of the model structure as shown in Fig. 9. The differ- 
ences of the electron temperatures can reach 100 K, the gas pressures can differ by 20% 
and the opacity averages, in particular the Rosseland mean, can differ by close to 60 %. 
These effects are much larger in the outskirts of the atmosphere than the changes in the 
structures caused by different TiO line lists (see Fig. 6). As a result, the use of the AMES- 
H2O line list also affects the optical spectra, causing weaker TiO bands than obtained with 
the MT-H2O line list. Models of early-type M dwarfs based on the AMES-H2O line list 
therefore systematically predict yet lower effective temperatures for a given star. 
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3.3. Combined effects 

In Figs. 10 and 11 we display a comparison between NextGen models (which use 
MT-H2O and SC AN-TiO) and models that use the AMES-H2O and AMES-TiO. The wave- 
length range of important filters and band identifications for TiO are given on the fig- 
ures. The TiO bands in the "AMES-atmosphere" are considerably weaker than those of 
the NextGen model spectrum as a result of the smaller oscillator strengths used, and 
the structural effects. On the other hand, the water bands are stronger in the AMES- 
atmosphere than in the NextGen model. This model has a relatively high temperature, 
thus the higher energy levels of the water molecule are relatively more important than for 
models with lower Teff (however, the concentration of water molecules is reduced so the 
overall water opacity is smaller than in the cooler models). 

To better judge of the impact of these opacity changes on the overall SED of M dwarfs 
in general, we have computed synthetic photometry as described in AH95 for three sets 
of models: (1) the NextGen grid based on J94 and MT-H2O (2) the AMES grid based on 
AMES-TiO and AMES-H2O opacities, and (3) the AMES-MT grid based on AMES-TiO 
and MT-H2O opacities. The results are compared to a photometric sample of M dwarfs 
(Leggett 1992) in Figs. 12, 13, and 14. Since M dwarfs form a tight sequence in optical 
VRI two-colors diagram despite the age and metallicity scatter of the sample (see AH95), 
this diagram imposes a strong constraint on model atmospheres. We find that models 
based on AMES-TiO opacities are systematically redder inV — R and V — I than models 
based on the J94 line list. The new models agree much better with observations and 
the new TiO data removes most of the discrepancy shown by the NextGen models in 
the lower main-sequence. Small remaining discrepancies may be attributed to the JOLA 
handling of VO and CaH which tend to overestimate slightly their opacities in the present 
models. Leinert et al. (1999) already studied the low resolution HST/FOS spectra of an 
M6 dwarf (LHS1070A) and found the AMES-MT models indeed agree quite well both 
with the observed SED and absolute fluxes within errors on the parallax of the system. 
They however noticed that some "continuum" flux excess remains important in the visual 
part of the SED (0.45 to 0.65 /xm). However, there is no a priori reason to assume that the 
a-f oscillator strengths are inaccurate, and these remaining problems could be related to 
other effects on the model structure. 

The use of AMES-H2O seems also to bring some improvements to the modeling of 
near-infrared colors. Fig. 13 shows that the late-type dwarfs can be better reproduced by 
the new water opacities than by the MT-H2O line list. However this diagram is sensitive 
to both gravity (lower gravity models loop lower) and metallicity, which makes it difficult 
to constrain the models on the adequacy of the water opacities used with them. Leggett, 
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AUard, & Hauschildt (1998) and Leggett et al. (1999) have already used the AMES models 
in their analysis of M dwarfs and brown dwarfs, and found an excellent general agree- 
ment of the predicted near-infrared SED with observations. However these analyses used 
the models to derive the parameters of the studied stars and brown dwarfs based on fits 
to the near-infrared SED or photometry and could not make an independent statement 
on the quality of the water line list. 

M dwarfs form again a sequence in the mixed-colors IJK diagram (Fig. 14), although 
less tightly than in the VRI diagram. Unresolved binary stars produce K-band flux ex- 
cess and lie below the sequence. H2 pressure-induced opacities depress the K-band flux 
of metal-depleted dwarfs so that they systematically lie above the sequence. But, as op- 
posed to the JHK diagram, this one is not particularly sensitive to gravity in M dwarfs 
which allows a sequence to be defined. Models should pass therefore through the bulk 
of early-type M dwarfs at } — K — 0.8, and follow a relatively / — K-insensitive sequence 
towards late-type dwarfs. We find that models based on the AMES-H2O opacities lie, as 
did the AH95 models before them, 0.2 magnitude in / — K to the blue of the observed 
sequence! And our tests show that this result is independent of the TiO opacities used. 
The NextGen models already reproduced perfectly the location of lower main-sequence 
stars in this diagram. And AMES-MT models computed using the AMES-TiO and MT- 
H2O line lists behave adequately both in the optical and infrared. Why? Perhaps the 
new water vapor line list is still not complete enough to high temperatures and lacks 
opacity in the /-bandpass i.e. around 1.3 /xm? Or would it have too much opacity in the 
K-window i.e. around 2.2 /xm? Until these questions can be answered, we hope that the 
two main grids of models we have computed (AMES and AMES-MT) will allow indepen- 
dent detailed confrontations to observations of cool stars that will locate more precisely 
the source of the problem (e.g. Leinert et al, in preparation). 



4. Summary and Conclusions 

A long standing problem with M dwarf models was that prior TiO line lists were 
incomplete to high temperatures. The use of "straight means" (AH95 models) helped 
by the coarseness of the treatment to block flux which otherwise escapes between lines 
in the incomplete list. But these models also blocked too much flux in most cases, and 
were only appropriate in late-type M dwarfs when TiO bands are already very strong. 
Clearly, a more complete line list was needed to model stars from the onset of TiO forma- 
tion to its gradual disappearance from the gas phase in brown dwarfs. The AMES-TiO list 
now serves beautifully this purpose. We find that the list provides more opacity in most 
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bands and suppress adequately flux between bands. The new, smaller oscillator strength 
values also play an important role in systematically assigning cooler models (at least for 
early type M dwarfs) to a given star, this way contributing to broader bands and lesser 
inter-band flux as well. These effects combine and should resolve most of the previously 
observed discrepancy between models and observations in the optical SED and photom- 
etry of M stars. Leinert et al. (1999) note however that flux excess remains substantial 
in the visual spectrum, suggesting some further incompleteness or fei inaccuracies of the 
new TiO in the a-f system. 

In order to better reproduce the observed (V-I) color indices, we had to retain in the 
present models the Davis, Littleton, & Phillips (1986) fei-values for the two reddest band 
systems: S and (p. For these two band systems, the theoretical estimates of Langhoff (1997) 
predicts an unexpectedly large fei-value ratio, while no laboratory estimate (Hedgecock, 
Naulin, & Costes 1995) are available to corroborate this. And we find, as did Alvarez & 
Plez (1998) in red giants for the 6 band, that models based upon the DLP86 fepvalues for 
these two bands reproduce adequately their observed depths in M dwarfs. 

The introduction of the AMES-H2O opacities bring solid improvements of the near- 
infrared SED of late-type dwarfs, but fails as the AH95 models did to reproduce ade- 
quately the J — K colors of hotter stars. Water vapor is a more important factor for the 
structure of the atmosphere than TiO because its overall opacity is larger and its lines are 
closer to the peak of the SED than the TiO bands, so the flux blocking effect of water vapor 
is more important for the temperature structure than that of TiO opacities for these low 
temperatures. Schwenke and collaborators at NASA AMES are preparing a new dipole 
moment function for H2O, which may change the high temperature high overtone water 
bands, and help resolve this discrepancy in the near future. 

Until a revised version of the AMES-H2O line list becomes available, we have there- 
fore generated two sets of model atmospheres for cool stars which allow to investigate 
these issues: the AMES grid based on the new TiO and H20 opacities, and the AMES-MT 
grid which rely on the AMES-TiO and MT-H2O opacities. 
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5. Tables 

Table 1: TiO fei values (fei = fvivi/qvivi) 



system 


lamO 


B90 


J94* 


DLP86 


L97 


AP98 


adopted 


alpha 


5170.7 


0.10 


0.17 


0.106 


0.105 


0.106 


0.105 


beta 


5605.2 


0.15 


0.28 


0.125 


0.176 


0.125 


0.176 


gamma' 


6192.5 


0.08 


0.14 


0.0935 


0.108 


0.0935 


0.108 


gamma 


7095.8 


0.09 


0.15 


0.0786 


0.092 


0.0786 


0.092 


epsilon 


8407.6 


0.0024 


0.014 


< 0.006 


0.002 


0.0023 


0.002 


delta 


8870.9 


0.02 


0.048 




0.096 


0.048 


0.048 


phi 


11044.8 


0.02 


0.052 




0.018 


0.0178 


0.052 



* laboratory values determined by HNC95 
B90: Brett (1990) 

DLP86: Davis, Littleton, & Phillips (1986) 
J94: Jorgensen (1994) 

HNC95: Hedgecock, Naulin, & Costes (1995) 
L97: Langhoff (1997) 
AP98: Alvarez & Plez (1998) 
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6. Figures 

Fig. 1. — Comparison between solar abundance M dwarf models calculated using AMES- 
TiO (full curves) and SCAN-TiO (dotted curves) in the blue spectral region. 
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Fig. 2. — Comparison between solar abundance M dwarf models calculated using AMES- 
TiO (full curves) and SCAN-TiO (dotted curves) in the red spectral region. 

Fig. 3. — Comparison between Teff = 2900 K; log(^) = 5.0 and solar abundance models 
calculated using AMES-TiO (full curves) and SCAN-TiO (dotted curves) based on the 
same set of TiO oscillator strengths (see "adopted" in Table 1). The models use otherwise 
identical opacities and parameters. Each model is fully converged, and the synthetic 
spectra are downgraded to a resolution of 20A. The positions of the TiO band heads are 
indicated according to table 8 of Langhoff (1997). 

Fig. 4. — Comparison between solar abundance M giant models calculated using AMES- 
TiO (full curves) and SCAN-TiO (dotted curves) in the blue spectral region. 

Fig. 5. — Comparison between solar abundance M giant models calculated using AMES- 
TiO (full curves) and SCAN-TiO (dotted curves) in the red spectral region. 

Fig. 6. — Comparison between solar abundance M dwarf models with Teff = 2500 K calcu- 
lated using AMES-TiO and SCAN-TiO. The differences are calculated in the sense AMES- 
TiO minus SCAN-TiO model. The bottom panel gives the results for the Planck-mean 
(full curve), J-mean (dotted curve). Flux-mean (dashed curve) and Rosseland mean dot- 
dashed curve) opacities. 

Fig. 7. — Comparison between solar abundance M giant models with Teff = 3000 K calcu- 
lated using AMES-TiO and SCAN-TiO. The differences are calculated in the sense AMES- 
TiO minus SCAN-TiO model. The bottom panel gives the results for the Planck-mean 
(full curve), J-mean (dotted curve). Flux-mean (dashed curve) and Rosseland mean dot- 
dashed curve) opacities. 

Fig. 8. — Comparison between solar abundance M dwarf models calculated using AMES- 
H2O (full curves) and MT-H2O (dotted curves) in the near-infrared spectral region. Both 
sets of models were calculated using AMES-TiO and iterated to convergence with their 
respective parameters. 
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Fig. 9. — Comparison between solar abundance M dwarf models with Teff = 2500 K 
calculated using AMES-H2O and MT-H2O. The differences are calculated in the sense 
AMES-H2O minus MT-H2O. The bottom panel gives the results for the Planck-mean 
(full curve), J-mean (dotted curve). Flux-mean (dashed curve) and Rosseland mean (dot- 
dashed curve) opacities. Both sets of models were calculated using AMES-TiO and iter- 
ated of convergence with their respective parameters. 

Fig. 10. — We compare the optical spectral distribution of a NextGen model with Teff = 
3500 K; log(g) = 5.0;[M/H] = 0.0 (dotted line) with a model converged on the same pa- 
rameters using the AMES-TiO and H2O line lists (full line). The positions of the TiO band 
heads are indicated according to table 8 of Langhoff (1997). The region of integration 
of standard optical broadbands are also shown for reference. The "AMES atmosphere" 
shows weaker TiO bands, principally due to the smaller oscillator strengths predicted by 
the Langhoff TiO model. 

Fig. 11. — Same as in Figure 10 for the near-infrared portion of the spectrum. The region 
of integration of the standard near-infrared broadbands are also shown for reference. The 
"AMES atmosphere" shows stronger H2O bands, especially in the trough of the bands, 
i.e., at 1.6 fxm and at 2.1 /im. But little or no changes are seen in the J bandpass region. 

Fig. 12. — The optical Cousins broad-band synthetic photometry of solar metallicity and 

fixed gravity (log(^) = 5.0) models of the NextGen grid (dotted line), AMES-MT grid 
(long-dash line), and AMES grid (full line) are compared to the photometric sample of 
Leggett (1992). This sample contains mostly M dwarfs and metal-depleted M subdwarfs 
of the solar neighborhood, and becomes scarce in the late-type dwarf regime. 

Fig. 13. — Same as Fig. 12 for near-infrared broad-band colors covering the water opacity 
range. Please note that the hot star tail of the sample, near H-K=0.1, is reproduced by the 
NextGen and AMES-MT models for the lower gravities predicted by evolution models 
for 5 to 10 Gyrs isochrones. Of the models shown, only the NextGen are grain-less, which 
explains their curling up at the low temperature end compared to AMES-MT models. 

Fig. 14. — Same as Fig. 12 and 13 for broad-band colors sampling side-to-side of the SED's 
flux peak. 



